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We develop and solve a continuum theory for the piezoelectric response of one dimensional nan- 
otubes and nanowires, and apply the theory to study electromechanical effects in BN nanotubes. We 
find that the polarization of a nanotube depends on its aspect ratio, and a dimensionless constant 
specifying the ratio of the strengths of the elastic and electrostatic interactions. The solutions of 
the model as these two parameters are varied are discussed. The theory is applied to estimate the 
electric potential induced along the length of a BN nanotube in response to a uniaxial stress. 
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Recently there has been interest in the physical prop- 
erties of materials that are obtained by wrapping a two 
dimensional lattice to form a one dimensional structure. 
When the lateral dimensions of the wrapped structure are 
reduced to the nanometer scale, the new periodicity in- 
troduced by this wrapping can have a profound impact on 
the physical properties of the system. The III-V boron- 
nitride (BN) exhibits a ground state polarization and a 
piezoelectric response (electric polarization linearly cou- 
pled to an elastic strain) when it is wrapped into a nan- 
otube, although the symmetry of the two dimensional BN 
sheet permits only a piezoelectric response. This system 
thus demonstrates the possibility of geometrical control 
of the piezoelectric response of a nanomaterial, provid- 
ing a new degree of freedom with significant potential for 
nanoscale electromechanical applications. 

The piezoelectric response of a macroscopic three di- 
mensional system is well described by the continuum 
Landau Devonshire model 0. However, this model docs 
not adequately describe the piezoelectric response of a 
one dimensional system because the ID Coulomb kernel 
is a short range potential (its Fourier transform diverges 
at small momentum q proportional to — log(g) rather 
than as l/q 2 as in three dimensions). Here we show that 
this short range kernel leads to new physical effects in ID 
systems such as nanotubes and nanowires. We find that 
the elastic response of a finite length piezoelectric nan- 
otube to an applied uniform mechanical load is in general 
spatially nonuniform. The spatial variation of the elec- 
tric polarization produces a bound charge density that 
penetrates from the tube end into the interior of the sys- 
tem. In addition, localized surface bound charges develop 
exactly at the tube ends. The partitioning of the bound 
charge into its "bulk" and "surface" contribution is then 
determined by a competition between the electrostatic 
and elastic interactions in the system. These effects are 
not described in any simple model for an electrostatic de- 
polarizing field for the tube, as would be possible for three 
dimensional samples in confined geometries (e.g. spheres, 
ellipsoids or thin films), and instead require an appropri- 
ate continuum theory to study the spatial variation of 
the electromechanical coupling. In this paper we develop 



this model and apply it to estimate the magnitude of ex- 
perimentally accessible electromechanical effects for the 
case of a BN nanotube. 

We consider the effect of an applied mechanical stress 
on a piezoelectric nanotube with radius R and length 
L. (The generalization to the case of a nanowirc with a 
solid core is straightforward.) A constant stress a(z) = f 
is applied to the tube and it induces a strain field rj{z) 
where z is the coordinate along the tube axis. This strain 
is linearly coupled to an electric polarization along the 
tube axis through a piezoelectric constant e s , so that 
P(z) = e s rj{z). In general an electric polarization of the 
nanotube can be produced by the uniaxial strain (exten- 
sion or compression along the tube axis) and by tube 
torsion. For example, for BN nanotubes the piezoelec- 
tric constants that couple to these two elastic strains are 
determined by the chiral angle that specifies the tube 
wrapping. For simplicity, in this paper we treat the case 
where only one elastic degree of freedom is coupled to 
the polarization. For the BN nanotubes this is realized 
by "zigzag" tubes that have a nonzero piezoelectric cou- 
pling only to uniaxial strain, or for "armchair" tubes that 
couple only through the tube torsion [2j ■ 

The free energy of the system can be written as a sum 
of its bare elastic and electrostatic contributions 

G — Gdastic + Gdcctrostatic — / r](z)o~(z) dz . (1) 
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Introducing the one dimensional clastic modulus C\ , the 
elastic energy is 

Gdastic = \C X [ V (z) 2 dz. (2) 

Here C\ = 2irRC2 where C2 is the elastic modulus 
for the two dimensional sheet (e.g. Cn for a uniaxial 
strain). C2 is the relevant intensive quantity (with units 
energy /area) that depends on the composition and struc- 
ture of the tube, and the one dimensional tube modulus 
C\ is an extensive quantity proportional to the tube ra- 
dius. 

Using the piezoelectric constant e s the electrostatic en- 
ergy can be expressed in terms of the spatial derivatives 
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of the strains as 
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Here e s is a one dimensional piezoelectric constant giving 
the strain-induced dipole moment per unit length (with 
units of charge) . For BN nanotubes of moderate radius it 
can be well approximated using the piezoelectric constant 
of an infinite two dimensional sheet e 2 mapped onto the 
tube circumference |2( ■ Thus for most physically relevant 
situations e s » 2nRe 2 - V(z — z') is a Coulomb kernel 
describing the electrostatic interaction between rings of 
charge centered at positions z and z' along the tube axis. 

It is useful to extract the dimensional dependence of 
these energies by expressing all lengths in units of the 
tube radius. Then by introducing the scaled variable 
£ = z/R, the aspect ratio Ai = L/ R, and another dimen- 
sionless parameter A2 = e\j{RC 2 ), we have 
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Here (3 = f / (211RC2) is the equilibrium strain in the 
absence of any electrostatic coupling. The competition 
between the electrostatic and elastic effects on the tube 
is therefore determined by the relative sizes of the scaled 
lengths Ai and A2. In the regime where A2 <C Ai the 
electrostatic interactions are negligible, and the equilib- 
rium strain is nearly constant along the length of the 
tube, albeit with a reduction near the tube ends. Con- 
versely when A2 Ai, the electrostatic interactions dom- 
inate and the system seeks to minimize the spatial deriva- 
tives of the strain, subject to the constraints of mechan- 
ical equilibrium. This leads to a situation where the 
bound charge density penetrates deep into the tube in- 
terior and the applied mechanical stress is nearly com- 
pletely balanced by the electrostatic interactions. Equa- 
tion J3J| also demonstrates that because of the scaling of 
the one-dimensional modulus with tube radius, the in- 
duced strains are inversely proportional to R, i.e. for 
a fixed aspect ratio large radius tubes are mechanically 
"stiffer". 

The equilibrium strain field is obtained by minimizing 
G and satisfies 

«=^M 2 r«fe (5) 



Possible stepwise discontinuities of rj(z) at the sample 
boundaries can be treated by explicitly writing 



(6) 



where r\i denotes the strain field in the sample interior. 
Then, by integrating Eq. JSJ by parts we get 
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In the following we will drop the subscript i, and the £' 
integrals are understood to extend over the interior of the 
tube. 

The solution of Eq. J7J requires a specification of an 
appropriate Coulomb kernel V(z). In the simplest con- 
tinuum model where V(z) describes the interaction of 
rings of charge, we have 
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This kernel diverges logarithmically at short range. Thus 
in this model a surface localized bound charge density 
p(0)5(z) at the tube end leads to a divergent electrostatic 
energy. Therefore the boundary condition for the strain 
field in this model becomes r/(0) = rj(L) — 0. This infinite 
energy is an artifact of the model and signals a breakdown 
of the continuum theory at short length scales. To avoid 
this difficulty, we instead use a "softened" kernel 
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which retains the correct long range behavior but remains 
finite for This softened kernel allows for the 

presence of a nonzero surface charge density e s rj(0)S(z) at 
the ends the tube. The value of a determines the energy 
"cost" of surface charges, and therefore controls the value 
of 77(0). Ultimately the distribution of the bound charge 
into its "surface" and "bulk" contributions is determined 
by a self consistent solution to the equation of state (J7J. 

Once the equilibrium charge density is determined the 
electric potential at position z along the tube is obtained 
by evaluating 
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After integrating by parts and transforming to scaled 
variables the potential is expressed 

u(o = 2^2 £ 9Vi l~ c \ (0 <e (11) 

Since the equilibrium condition in Eq. (JJJ requires that 
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the potential can be expressed as an integral over the 
strain field 77, and we obtain 
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(13) 
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where we have set the electric potential to zero at the 
center of the tube. 

The equilibrium solutions for the strain field rj{z) (Eq. 
EJ, the bound charge density p — ~e s drj/dz and the elec- 
trostatic potential U(z) (Eq. [T3*j l are displayed in Fig. |QJ 
for the model in the regime where the elastic interactions 
dominate (Aa <C Ai) and in Fig. (j2J) for the regime where 
the electrostatic interactions dominate (Ai <C A2). The 
softened kernel was used with a = .001. The data in 
these plots was obtained by finite clement analysis. The 
tube length was divided into 1000 unequally spaced bins 
(1/3 evenly spaced in the first 10% of the tube, 1/3 evenly 
spaced in the last 10%, and 1/3 evenly spaced across the 
middle 80% of the tube). Bin number was increased in 
steps of 100 until successive solutions differed by less than 
1.2%. The slowest convergence occurs at the end of the 
tubes, with the solution differing by less than .2% on the 
last iteration over the middle 99.8% of the tube. 

In the elastically dominated limit, the equilibrium 
strain is nearly constant in the interior region of the tube 
where 77 « (3. The electrostatic interactions effectively in- 
crease the elastic moduli near the surface region, leading 
to a suppression of the strain. This suppression implies a 
spatial variation of the polarization and a volume bound 
charge density that is obtained self consistently by solv- 
ing Eq. (0). Finally, we see that the equilibrium strain 
is in general nonzero at the tube ends. Thus, in addition 
to the continuous volume bound charge density, a true 
surface localized bound charge is obtained at the termi- 
nations. In a long (semi-infinite) tube the equilibrium 
strain asymptotically relaxes to its elastic limit following 
a power law, e.g. 1 — r)(£)/P oc (7/£) 2 where 7 = y/2n\2. 
The distribution of the bound charge is plotted in the 
central panel, which shows that the main effect of the 
electrostatic interaction is to "spread" the bound charge 
density near the tube ends. In the elastically dominated 
limit most of the potential drop across the tube, shown 
in Fig. Q-d), occurs in the regions near the tube ends. 

This behavior is contrasted with the situation for the 
electrostatically dominated limit shown in Fig. @. In 
this regime, the system seeks to reduce its strain gradi- 
ents subject to the constraint of mechanical equilibrium. 
This leads to a state in which the bound charge extends 
along the entire length of the tube. The numerical solu- 
tion for the equilibrium strain shown in Fig. J2J-a) is rea- 
sonably well approximated using only the lowest allowed 
Fourier component with the ends of the tube clamped. 
The overall magnitude of the strain is also strongly re- 
duced, since it is limited by the electrostatic coupling 
strength instead of the bare elastic constant. The bound 
charge density vanishes at the center of the tube by sym- 
metry and varies approximately linearly as a function of 
position, with logarithmic corrections close to the ends of 
the tube. Interestingly, Eq. (fT2"|) requires that the elec- 
tric potential is a linear function of position in the limit 
that i](z) <§C /3, as is shown in Fig. 10-d). 
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FIG. 1: Equilibrium solutions for a piezoelectric tube in 
the elastically dominated limit, (a) Equilibrium strain profile 
from Eq. 0. The strain is well approximated by the elastic 
limit over most of the tube, (b) Inset, strain value at tube 
boundary. Notice r/(0) 7^ with our choice of V. (c) The 
bound charge density, p = —e s dr\jdz. The arrows represent 
delta functions in the charge density, (d) Equilibrium poten- 
tial from Eq. 1131 . The graphs were generated using (3 — 1, 
Ai = 1000, A2 = 10/-7T, and e s = 1 (charge). 



Since the one dimensional Coulomb kernel is a short 
range interaction it is useful to consider the limit where 
the kernel is a contact interaction V(£ — £') = S(£, — £'). 
Then the free energy in Eq. I0J reads 



G/{ttR 2 C 2 ) = 



(77(0) 2 " 2/377(0 + 27TA; 



(14) 

Up to an additive constant, this describes the free energy 
of a superconducting disk, with penetration depth \/2n\2 
immersed in a perpendicular field Hq — —f3. The equi- 
librium solution to this model is obtained analytically as 



v(0 = 01- 
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FIG. 2: The equilibrium solutions for a piezoelectric tube in 
the electrostatically dominated limit, (a) Equilibrium strain 
profile from Eq. |0. Notice that the maximum value of the 
strain is about 1/5 of the elastic limit, (b) Inset, strain value 
at tube boundary. Notice 77(0) 7^ with our choice of V. (c) 
The bound charge density, p = —e s dn/dz. The arrows rep- 
resent delta functions in the charge density, (d) Equilibrium 
potential from Eq. 1131 . The graphs were generated using 
= 1, Ai = 1000, A 2 = 10000/tt, and e s = 32 (charge). 



Thus the strain profile in the elastically dominated limit 
shown in Fig. (JTJa) give the piezoelectric analog to the 
Meissner expulsion of an applied field (stress) and the 
electrostatic regime of Fig. (PJ-a) corresponds to the limit 
of a small superfiuid density where the penetration depth 
is of order the length of the system. For the contact inter- 
action model the induced potential is simply proportional 
to the bound charge density, so that a measurement of 
the spatial variation of the potential probes the piezo- 
electric analog of the London penetration depth. Note 
however that the physical kernel retains an asymptotic 
l/£ tail leading to a power law decay of the strain field. 
Recently, BN nanotubes were predicted to exhibit a 



significant piezoelectic response 0, Q . Here we predict 
the piezoelectric response of a "zigzag" BN nanotube un- 
der uniform longitudinal tensile or compressive stress. 
The 2D elastic constants for various chiral BN tubes were 
calculated in Ref. |3 here we use an average value of 
C2 = -30 TPa nm . Ref. 2 gives the piezoelectric con- 
stant for a BN sheet as e-2 — .12 e/Bohr. Experimentally 
accessible forces are on the nano-Newton scale; we arbi- 
trarily assume a force of 1 nN. Typical tube dimensions 
are R = 1 nm and L = 1 jum. 

Using these values we find that Ai = 1000, A 2 = 
4.4 x 10~ 3 , and f3 — 5.3 x 10~ 4 . These values give an 
asymptotic decay constant of 7 = .17, and the strain to 
relaxes to its elastic limit very quickly. Our numerical 
solution shows the equilibrium strain is within 1% of the 
elastic limit at a depth of 3i?. Thus the bound charge 
density, although continuous, barely penetrates the bulk 
of the tube, and the entire potential drop occurs near 
the ends. Applying the continuum theory to this sys- 
tem and converting to SI units we obtain an estimate of 
.48 V across the tube for an applied tension / = 1 nN. 
It would be interesting to experimentally measure this 
strain induced electrostatic potential and to resolve its 
spatial variation near the tube ends using scanning probe 
potentiometry. 

The theory for a solid core nanowire follows the same 
formalism, but unlike the nanotube, in a nanowire the 
corresponding A2 is fixed entirely by the material param- 
eters, independent of R. This occurs for a nanowire be- 
cause the extensive variables, e s and Ci , are proportional 
to the cross-sectional area rather than the circumference. 
They are related to the relevant intensive variables by 
e s = uR 2 e^ (where is the 3D piezoelectric constant 
in units of charge/area) and C\ = ttR 2 Cs (where C3 is 
the 3D elastic constant in units of pressure). The equi- 
librium strain then satisfies Eq. Q with A2 = e\j(2C^) 
and (3 = f / {ttR 2 C^). Thus two nanowires with the same 
aspect ratio but different radii will have the same scaled 
strain, 77//?, unlike hollow nanotubes where larger radius 
tubes have smaller A2 and thus are elastically dominated. 
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